Vortices and turbulence in trapped atomic condensates 



cn 

o 

u 

< 
0\ 



c/3 

I 



I 

o 
o 



> 

m 
en 
in 

^' 

o 
en 

> 

X 



Angela C. White/ '[j Brian P. Anderson,^ and Vanderlei S. Bagnato'^ 

'Joint Quantum Centre (JQC), Durham- Newcastle, School of Mathematics and Statistics, 

Newcastle University, Newcastle upon Tyne, NEl 7RU, UK. 

^College of Optical Sciences, University of Arizona, Tucson, AZ 85721, USA 

^ Instituto de Fisica de Sao Carlos, Universidade de Sao Paulo, 

Caixa Postal 369, 13560-970 Sao Carlos, SP, Brazil. 

After over a decade of experiments generating and studying the physics of quantized vortices 
in atomic gas Bose-Einstein condensates, research is beginning to focus on the roles of vortices in 
quantum turbulence, as well as other measures of quantum turbulence in atomic condensates. Such 
research directions have the potential to uncover new insights into quantum turbulence, vortices and 
superfluidity, and also explore the similarities and differences between quantum and classical turbu- 
lence in entirely new settings. Here we present a critical assessment of theoretical and experimental 
studies in this emerging field of quantum turbulence in atomic condensates. 



INTRODUCTION 

Since Onsager's groundbreaking theoretical work link- 
ing turbulence and point vortex dynamics in a two- 
dimensional (2D) fluid [T], it has been hoped that the 
simple nature of quantum vortices in superfluids will aid 
in understanding the nature of turbulence. After many 
years of research with supcrfluid helium systems, the field 
of quantum turbulence (QT) is now well established, and 
has led to numerous new insights and developments re- 
garding QT and the universality of turbulence [2]. The 
discovery of links between classical turbulence and QT 
remains a strong motivating factor for QT research, par- 
ticularly in the emerging field of QT studies with Bose- 
Einstein condensates (BECs) . BECs present a new plat- 
form for QT studies due to their compressibility, weak 
interatomic interactions, and availability of new experi- 
mental methods for probing and studying superfluid flow 
[3]. The relationship between QT and vortex dynamics 
in these systems is consequently an inherently interesting 
new research topic as well. 

Classical turbulence is composed of eddies of contin- 
uous vorticity and size and it is necessary to solve the 
Navier-Stokes equation to mathematically describe vis- 
cous fluid dynamics [4j. For turbulent fluid flow, which 
consists of scale-invariant flow dynamics across a wide 
range of length scales, this procedure becomes difficult to 
tackle from first principles. In comparison, QT is com- 
prised of vortices of less complexity, each with a localized 
and well-defined vortex core structure and quantized cir- 
culation. Superfluid flow is inviscid and vortices can- 
not decay by viscous diffusion of vorticity: a quantized 
vortex cannot simply "spin down" and dissipate energy 
via viscosity in the same way a classical vortex can. In- 
compressible kinetic energy is instead diffused through 
emission of sound waves and then dissipated due to the 
presence of a thermal cloud in BECs or the normal fluid 
component in superfluid He. 

Despite the differences arising from the nature of vor- 
tices, classical and quantum turbulence share profound 



similarities that underscore the universality of turbu- 
lence. We briefly illustrate this idea with the structure of 
kinetic energy spectra in three-dimensional (3D) turbu- 
lence. At locations in the fluid far from vortex cores, 
and for length scales greater than the average inter- 
vortex spacing, vortex core structure is unimportant and 
quantized vortex lines are analogous to vortex fllament 
lines of an Euler fluid. Under these conditions, classical 
and quantum turbulence are known to possess similar 
macroscopic and statistical properties. The most strik- 
ing similarity is the existence of the same Kolmogorov 
spectrum in the inertial range for 3D turbulence. For 
non-equilibrium steady-state forced turbulence, at length 
scales larger than the distance given by the average inter- 
vortex separation and smaller than the scale correspond- 
ing to energy injection, the incompressible kinetic energy 
spectrum scales with wavenumber k as E{k) ^ k~^'^ 
[5]. This scaling is thought to arise from a Richardson 
cascade process as large vortices break up into smaller 
and smaller vortices until at very small scales energy is 
dissipated. The same classical Kolmogorov scaling has 
been verifled numerically and experimentally in turbu- 
lent superfluid '*He and ^He-B [51-fTT] and has also been 
established numerically as a feature of trapped [121 ES] 
and homogeneous [111 [15] atomic condensates. There is 
thus reason to believe that a quantitative understanding 
of aspects of turbulence in one system, even a quantum 
fluid, may aid in the general understanding of the sub- 
ject. 

Among the central issues in the development of an un- 
derstanding of QT is therefore understanding how quan- 
tized vortices move about in a plane in 2D QT, or how 
they bend, interact, and create complex tangles in 3D 
QT. The relationship to vortex dynamics and statistical 
measures of turbulence, such as energy and velocity spec- 
tra, is also an ongoing and active research topic. While 
the physics of such processes and the links between the 
various aspects of the problems may remain intricate, 
the reduction of vorticity to well-defined cores enables 
new approaches to studying QT generally unattainable 



in classical turbulence. This article is aimed at present- 
ing a brief overview of theoretical, numerical, and experi- 
mental progress — as well as potential capabilities — for 
achieving a deeper understanding of QT by tackling these 
central issues. We present atomic BECs as a new and 
distinct tool that will contribute towards an extended 
understanding of QT and vortex dynamics through both 
theory and experiment. Firstly we discuss the unique 
aspects of turbulence in BECs in comparison to super- 
fluid ^He and •^He-B. We examine the recent theoretical 
and experimental advances in the BEC turbulence field 
and detail the new regimes of turbulence made accessible 
in BECs. Finally we overview experimental measures of 
turbulence and future directions for the field, focusing on 
progress in both 3D and 2D QT. 



ATOMIC BECS: BEYOND SUPERFLUID 
HELIUM 

Studies of quantum turbulence have historically been 
restricted to superfluid "^He and ^He-B. These systems 
have a huge range of accessible length scales and as a 
result, turbulent vortex tangles can consist of hundreds 
of thousands of vortices. The vortex tangles are well 
separated with typical inter-vortex separation distances 
of I ^ 10~'*m. Vortices have small vortex core diame- 
ters, ^ ^ 10~^°m for superfluid ''He, and consequently 
the turbulent tangles in '*He and ■^He-B are characterised 
by a large ratio of inter-vortex spacing to vortex core 
radius, on order of 10^ — 10^. The intrinsic superfluid 
parameters such as atom-atom interaction strength are 
fixed, and the superfluids are homogeneous and of con- 
stant density (i.e., incompressible). Controlling single- 
vortex dynamics in a turbulent ^He and '^He-B superfluid 
is exceptionally challenging, and probing the behaviour 
of turbulence at small scales remains an open problem 
from an experimental standpoint. Finally, due to the 
strong interactions between atoms in liquid and super- 
fluid helium, simulations of wave function dynamics are 
generally only qualitatively accurate. 

In contrast, trapped atomic condensates are rarely ho- 
mogeneous and have a smaller range of accessible length 
scales over which vortex dynamics can be probed. With 
the exception of shallow traps and nominally hard-wall 
conflning potentials where the condensate is to a good 
approximation homogeneous, for realistic trapped sys- 
tems the condensate density is non-uniform, with the 
non-uniformity arising as a consequence of the form of 
the trapping potential and the compressibility of BEC 
systems. It is this compressibility that also lends signif- 
icant new approaches to the study of QT. Also because 
of this, manipulating the trapping potential allows one 
to manipulate the condensate homogeneity throughout 
the trapping region. The typical vortex core diameter 
in atomic BECs is much larger than that in superfluid 



He, being on order of the coherence or healing length, 
^ 0.5/xm. Turbulent BECs consist of small numbers of 
vortices (generally no more than ^ 100) that are less 
sparsely separated (ratio of inter- vortex spacing to vor- 
tex core radius on the order of 10). From a theoretical 
standpoint, the models describing vortex tangles in su- 
perfluids are amenable to numerical simulation, and for 
BECs are well established and quantitatively accurate at 
both zero and non-zero temperatures [THl HZ] . 

BECs also permit adjustment of intrinsic atomic prop- 
erties that lead to macroscopic flexibility not readily 
achieved in other superfluid systems. The strength of 
atom-atom interactions can be controlled and even driven 
from attractive to repulsive by tuning an external mag- 
netic field around a Feshbach resonance [18]. Individ- 
ual vortex dynamics and position can be well controlled 
[191 120] . and combined with a wide variety of imaging 
techniques opens up the field to investigations of few- 
vortex systems and the relationship of chaotic vortex 
dynamics to turbulent flow. The tunable features of 
trapped atomic condensates further extends to their di- 
mensionality. Highly oblate condensates, with vortex dy- 
namics well into the 2D regime can be routinely created. 
This makes BECs the first system for which 2D QT is 
experimentally accessible [20l - f23] . We note that for such 
studies, quasi-2D BECs that have chemical potential less 
than the most strongly confining dimension's mode spac- 
ing are not needed. Rather, highly oblate but never- 
theless 3D BECs may suppress superfluid flow along the 
tight-confining direction enough that 2D superfluid vor- 
tex dynamics and 2D turbulence are obtained [21]. We 
will henceforth refer to such BECs as 2D, and the asso- 
ciated conditions, where vortex dynamics occur within a 
plane, as being 2D. In addition to creating conditions for 
studying both 2D turbulence and 3D turbulence, BEC 
trapping parameters are extremely flexible, opening up 
the possibility of investigating transitions between 2D 
and 3D QT. 



Theoretical formalism 

The dynamics of quantized vortices in zero- 
temperature BECs with scalar order parameter are 
described by evolution of the Gross-Pitaevskii equation, 
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V^V (r) +F (r) V (r) + UqIiJj (r) fi/j (r) 



(1) 

Here i/" is a classical mean-field wave function repre- 
senting the condensate trapped by a potential V (r) . 
Uq = ATTfi?a/m describes interactions between bosons of 
mass m in the condensate, where a is the s-wave scat- 
tering length, and h is the reduced Plank's constant [3]. 
For simplicity, we consider BECs in cylindrically sym- 
metric harmonic trapping potentials. When the axial 



trapping frequency is much greater than the radial trap- 
ping frequency, w^ ^ a;,-, the condensate wiU be highly 
oblate, enabling studies of 2D QT and vortex dynamics 
limited to the radial plane. For lj^ ^ i^r, 3D QT may 
be explored. The dynamics of condensates at non-zero 
temperatures can be accurately described by applying 
classical field methods [17 . For a description of non-zero- 
temperature models that could be applied to superfluid 
turbulence, we refer the interested reader to the reviews 

[Him- 



Compressible and Incompressible kinetic energy 

Decomposing the condensate kinetic energy into com- 
pressible and incompressible parts is a useful technique 
frequently applied to analyse how kinetic energy due to 
vortex lines and sound is distributed over length scales 
[5]. This is performed by defining a density weighted ve- 
locity field T = v^v, in terms of the superfluid velocity, 
V = — V0, where n and are the position-dependent 
condensate density and phase profiles respectively. Uti- 
lizing the fundamental theorem of vector calculus, we can 
write T = T* -f Y"^, where T* is an incompressible (i.e., 
divergence-free) component satisfying V • T' = 0, and T'^ 
is a compressible component for which V x T^ = 0. The 
compressible and incompressible kinetic energy spectra 
E^kinik) are defined by 
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(2) 



fc| is a wavenumber and El^^ and E^^^ are the 



where k 

total incompressible and compressible kinetic energies in 
the system (per unit mass). The flow of incompressible 
kinetic energy across wavenumbers and its relationship to 
vortex dynamics is one of the central and critically im- 
portant and challenging issues in classical turbulence as 
well as QT. An understanding of spectra in 2D and 3D 
QT thus relates directly to the dynamics and distribu- 
tion of quantized vortices. Considerable effort has been 
recently concentrated on constructing the incompressible 
kinetic energy spectrum and also the angle-averaged mo- 
mentum and velocity distribution in order to determine 
their scaling in relation to the distribution of vortices in 
2D [23H51] and 3D QT in condensates ^. 



Quantum Pressure Energy 

Contributions to the quantum pressure energy, 
1 



E, 



IP- 2 I dr|V\Air , 



(3) 



arise only where the condensate density varies sharply, 
such as at vortex cores. This quantity may provide a 



useful theoretical measure of the total vortex number in 
2D condensates. For 3D condensates, quantum pressure 
energy scales with the vortex line length. While evaluat- 
ing incompressible, compressible and quantum pressure 
energies and spectra is useful from a theoretical point 
of view, it remains an open problem how to measure 
these individual components experimentally. However, 
velocity correlations and the angle-averaged momentum 
spectrum are quantities potentially experimentally acces- 
sible in atomic BECs and we outline some measurement 
prospects later in the article. 



VORTEX GENERATION 

A plethora of experimental methods are available to 
induce vortices in BECs |30j. We highlight a non- 
exhaustive selection of techniques that allow the prepa- 
ration of well defined initial states from which to inves- 
tigate the evolution or decay of turbulence. Determin- 
istic induction of vortices into the condensate at pre- 
cisely defined positions can be achieved by the controlled 
methods of phase imprinting, where the phase profile of 
the condensate is engineered [211 [31] , or by the coherent 
transfer of orbital angular momentum to a condensate 
by a two-photon stimulated raman process [55H55] . The 
creation of non-equilibrium vortex states with arbitrary 
winding has been beautifully demonstrated by the trans- 
fer of orbital angular momentum from a holographically 
produced light beam [37]. Applying this technique has 
the potential to create arbitrarily complex initial vortex 
distributions and even vortex knots in BECs [3H1 135] . 
Alternatively, vortex knots could be imprinted in BECs 
applying the techniques demonstrated in classical fluids, 
through the acceleration of shaped hydrofoils [40] , made 
by a 3D laser structure or a shaped nano-tube in a con- 
densate. 

Laser stirring can also be applied to create vortices in 
atomic BECs in a deterministic manner [U]. The dis- 
tribution of vortices throughout a condensate is known 
to depend on the path of the laser stirrer [42]. For 
a 2D condensate, this means the path the laser stirrer 
takes through the condensate can create clusters of like- 
signed vortices as shown in Fig. 1(A), or a more ran- 
dom distribution of vortices of differing sign. For a 3D 
condensate, this implies that a stirring path might be 
optimised to generate well-distributed vortex configura- 
tions or more polarised tangles, where the vortex dis- 
tribution is aligned preferentially along a particular di- 
rection. Reeves et al. [13] have shown that the potential 
strength of a laser stirrer and its speed can also be chosen 
such that the laser stirrer sheds single dipoles, clusters, 
or oblique solitons in a trapped 2D BEC. These ideas 
should also extend to stirring or flow past an obstacle in 
a 3D condensate. 

In addition to laser stirring, combined rotation and 



precession around three cartesian axes has been shown 
to create isotropic vortex tangles in 3D condensates [13] . 
Reducing the rotation and precession of the condensate 
to only two directions decreases the degree of isotropy 
of the resulting vortex tangle and could also be applied 
to create polarised distributions of vortices [H]. The 
methods of laser stirring, phase imprinting, transfer of 
orbital angular momentum and combined rotation and 
precession which have been highlighted here, have the 
additional advantage of producing only small amounts of 
phononic excitations, where the acoustic energy density 
is much less than the incompressible energy density; see 
[43] for further discussion of this relationship. 



CHAOS AND FEW- VORTEX DYNAMICS 

The ability to create well-defined initial distributions 
of vortices [23 opens up the possibility of directly prob- 
ing systems of few vortices and their resulting dynamics. 
Previous experiments have observed the precession of sin- 
gle filled vortices in trapped condensates [U], and the 
motion of vortex dipoles deterministically created prior 
to expansion of a sequence of BECs [31]. Recent ex- 
perimental advances have demonstrated measurements of 
few- vortex dynamics within a single BEC [45l |46] . This 
method is accomplished by allowing only a small frac- 
tion of the atoms in a BEC to be imaged after expan- 
sion, and sequential images from multiple expansion steps 
allows determination of vortex dynamics. This tech- 
nique could be incorporated into future measurements 
of chaotic vortex dynamics. With the future develop- 
ment of real-time in situ imaging of vortex dynamics it 
may become possible to directly observe vortex dynam- 
ics, including vortex-vortex annihilation and reconnec- 
tion events directly within BECs. The motion of four 
point vortices in a plane can be chaotic [47ti49] and sim- 
ilar chaotic dynamics are expected to be observable in 
oblate trapped BECs. The number of vortices that de- 
termines the crossover from chaotic to turbulent vortex 
dynamics and the role of chaos in turbulence [50] remain 
open questions that atomic BECs may be able to address. 



THREE-DIMENSIONAL TURBULENCE 
Non-equilibrium steady-state turbulence 

One of the defining features of continuously forced tur- 
bulence in a bulk 3D fluid is the existence of a direct 
Kolmogorov cascade corresponding to the conserved flow 
of energy from the forcing scale to smaller length scales. 
Although the direct Kolmogorov energy cascade has been 
established numerically for trapped atomic condensates 
[m [T3] , it is yet to be confirmed experimentally. Indeed 
for a realistic tangle of vortices, where the BEC diameter 



is typically less than ~ 100/im, the small range of length 
scales present may prohibit its detection. On the other 
hand, this implies trapped condensates may be a suitable 
system in which to determine the lower bound on the sys- 
tem size and vortex line density for which Kolmogorov 
scaling may be observed. Although the k^^/^ scaling 
of the incompressible energy spectrum is established in 
superfluid He systems and is thought to arise due to a 
Richardson cascade process [51-fTUl [5T] , a corresponding 
break up of large eddies into smaller and smaller eddies 
has not been theoretically or experimentally verified. If 
a Richardson cascade process is present in QT, it may 
be directly observable in trapped atomic condensates as 
detection schemes that experimentally image vortices di- 
rectly are feasible. We outline such techniques near the 
end of this article. Tracking the length of individual vor- 
tices is also possible numerically, and may help answer 
questions about the relationship of an energy cascade to 
a Richardson cascade. Establishing the existence or ab- 
sence of a Richardson cascade process coupled to a direct 
Kolmogorov cascade is perhaps one of the most exciting 
prospects in future experimental and numerical studies 
of turbulent vortex tangles in trapped atomic BECs. 



Small-scale dynamics and dissipation 

On length scales smaller than the typical inter-vortex 
spacing, quantum turbulence is quite different from clas- 
sical turbulence. Characteristics that are dependent on 
the vortex core structure and microscopic nature of the 
fluid are expected to be unique for quantum fluids. One 
such property is the velocity statistics of turbulent quan- 
tum fluids. The velocity components of trapped and 
homogeneous turbulent BECs in 2D and 3D, as well as 
turbulent superfluid He, follow power-law like behaviour 
[HI [52H51] , in stark contrast to gaussian velocity statis- 
tics of turbulent classical fluids [55H57] . This difference 
in the velocity statistics of quantum fluids arises from 
the singular nature of quantized vorticity and the l/r 
velocity held of a quantum vortex [531 ISH] • Measurement 
techniques in superfluid *He can also probe the 'quasi- 
classical' limit, that is they do not measure length-scales 
smaller than the average vortex separation distance. In 
this limit, velocity statistics lead to a gaussian distribu- 
tion [inilS^. This is another demonstration of features of 
classical turbulence emerging for quantum turbulent sys- 
tems, at large scales, where the microscopic vortex core 
structure is not probed. 

In 3D turbulence, the direct Kolmogorov cascade is 
coupled to dissipation of incompressible kinetic energy at 
small scales. In superfluid He, at zero temperature this 
sink of incompressible kinetic energy is facilitated by a 
Kelvin wave cascade process at scales smaller than the av- 
erage inter-vortex separation distance. The Kelvin-wave 
cascade process occurs when the helical perturbations on 



vortex lines cascade to shorter and shorter wavelengths, 
until at small scales phonons are excited, dissipating the 
kinetic energy of vortices [60j . For zero temperature con- 
densates, such a Kelvin wave cascade process could also 
provide a sink of incompressible kinetic energy, in par- 
ticular as Kelvin waves can be generated by vortex re- 
connection events (611 that are more frequent due to the 
smaller inter- vortex spacing in 3D condensates. As the 
collision of finite amplitude sound waves, or rarefaction 
waves with vortex lines can also generate Kelvin waves 
[62j . they may feature prominently in tangles generated 
by methods that induce a significant quantity of sound. 
Whether the Kelvin wave cascade process is driven by 
local or non-local transfer of energy remains an open and 
contentious question [55^70] . However, Gross-Pitaevskii 
equation simulations for Kelvin waves in a homogeneous 
system find scaling consistent with that predicted for a 
non-local energy cascade caused by 4-wave interactions 
[71] . Numerical simulations of few- vortex systems in 3D 
trapped condensates at non-zero temperature have shown 
Kelvin waves to be important in facilitating the decay of 
vortices by their movement out of the trap j24| . This 
work also showed that Kelvin mode excitations can be 
effectively frozen out by altering the condensate dimen- 
sionality through flattening the condensate, or increasing 
its confinement along the vortex line direction |24| . 

In addition to the dissipation of incompressible kinetic 
energy through a Kelvin wave cascade process, for finite 
size trapped condensates, the dissipation of incompress- 
ible kinetic energy by sound generated through reconnec- 
tion events [22j [73] will also be important. It has been 
suggested that the break up of vortex loops into smaller 
vortex rings could form a self-generating reconnection 
process to smaller and smaller scales until they are dis- 
sipated by self-annihilation and sound waves [TJJ [75] . 
Whether sound produced by a Kelvin wave cascade pro- 
cess, or from reconnection events, is the dominant process 
dissipating incompressible kinetic energy in zero temper- 
ature trapped atomic condensates, and how this scales 
with condensate size, vortex line density, and the polarity 
of vortex tangles, remain open questions to be addressed. 




FIG. 1: A). Density profile of a 2D condensate where a clus- 
ter of like-signed vortices has been induced by stirring with 
a laser paddle as in |42]. Vortices of positive and negative 
winding are denoted by magenta -I- and cyan — symbols re- 
spectively. B). 3D decaying turbulence: distribution of vortex 
lines (magenta) introduced by phase imprinting a lattice of 17 
straight vortices in a harmonically trapped condensate [53) . 
The condensate edge defined by a density that is 25% of the 
peak density is depicted in cyan. 



[55] . No distinction between scalings of vortex line-length 
decay predicted from 'Vinen' or 'ultra-quantum' mod- 
els of turbulence or 'quasi-classical' models of turbulence 
describing decaying turbulence in superfluid He was ob- 
served. 'Ultra-quantum' models of turbulence describe 
random distributions of vortices in superfluid He tur- 
bulence, with vortex decay facilitated by a Kelvin-wave 
cascade process corresponding to a decay of vortex line- 
density L a.s L ^ t~^ [751177]. For turbulence in super- 
fluid He where a Kolmogorov cascade features, vortices 
are locally aligned [75], and the decay of line density is 
observed to scale as L ~ i"^/^ [73 [771 [75]. A unique 
scaling of the decay of vortex line-length, dictated by the 
length-scale at which turbulence was forced, may feature 
in the decay of turbulent tangles in trapped condensates. 
As the number of vortices in a turbulent BEG is much 
less than in turbulent superfluid He, yet vortex tangles 
are typically more dense, the processes dictating vortex 
decay are expected to be strongly influenced by the den- 
sity of vortices and trapping inhomogeneity. 



Decay of turbulence 



In addition to forced 3D turbulence, it is also possible 
to investigate decaying QT in trapped 3D BECs. The 
decay of turbulence is expected to be a complex pro- 
cess facilitated by vortex-vortex reconnections with the 
rate of decay of vortex line density also influenced by 
the condensate homogeneity and temperature. A numer- 
ical study of the decay of a vortex tangle induced from 
straight vortices imprinted along each cartesian direction 
employed a phenomenological model of temperature and 
found vortex-line length to decay faster with increased 
dissipation, qualitatively modelling higher temperatures 



Experiments on three-dimensional QT 

Between the years 1999 and 2009, numerous aspects of 
the physics and characteristics of quantized vortices were 
studied in experiments with superfluid atomic gases [3D] . 
These investigations included methods of vortex gener- 
ation, imaging, and manipulation; properties of highly 
rotating condensates; vortex interactions, structure, sta- 
bility; and relationships to superfluidity. However, few 
references to turbulence can be found in this body of 
work, and appear primarily as brief qualitative state- 



ments given to refer to disordered distributions of vor- 
tices. Such observations were made in the contexts of 
rotating single-component [80l [81] and two-component 
[52] BECs, and with fast laser sweeps through a BEC 

m- 

In arguably the first experiment principally aimed at 
studying a turbulent tangle of vortices in a 3D BEC, 
a non-uniform distribution of vortices was created in a 
cigar-shaped condensate by exciting a surface mode in- 
stability induced by an applied external oscillatory po- 
tential [m [85]. Varying the strength and duration of 
the perturbing potential increased the number of vor- 
tices nucleated, enabling the observation of few-vortex 
configurations [86J through to a tangled distribution of 
vortices. In the regime of many vortices, it is difficult to 
identify individual vortex cores, particularly with detec- 
tion of an expanded BECs column density distribution. 
In these experiments on 3D QT, the turbulence was in a 
regime with few enough vortices that absorption images 
of expanded condensates showed the presence of vortex 
lines of varying orientations. A tangled vortex configu- 
ration was then inferred from the absorption images. On 
release of the trap, the condensate aspect ratio (defined 
by the measured ratio of BEC width to height) was con- 
served and self-similar expansion was observed [H] . This 
behaviour is in stark contrast to the typical inversion of 
the aspect ratio observed on expansion of a vortex-free 
cigar shaped condensate and is thought to be an indica- 
tion of a tangle of vorticity throughout the condensate. 
This conclusion is supported by theoretical investigations 
that studied the influence on the expansion of a cigar- 
shaped cloud arising from the presence of randomly dis- 
tributed vorticity [57] and distributions of vorticity with 
preferred directions in a hydrodynamic framework [55] . 
However, full numerical simulations going beyond the hy- 
drodynamic limit and including quantum effects are still 
required. Measuring the condensate momentum distri- 
bution could potentially also provide useful information 
about the vortex distribution throughout the condensate 
[25j . Measurements towards this aim are underway at 
the Unversidade de Sao Paulo. 



TWO-DIMENSIONAL TURBULENCE 

The existence of the direct Kolmogorov cascade for 3D 
superfluid turbulence suggests that large scale features 
of 2D turbulent classical flow [89H9T] may also crossover 
to 2D turbulence in quantum fluids. One of the charac- 
teristic features of classical 2D turbulence is the possible 
appearance of an inverse energy cascade (lEC) in which 
energy flow across length scales may occur in a direction 
opposite that of 3D turbulence; that is, with kinetic en- 
ergy injected at small length scales, energy flows towards 
larger length scales in an inertial range free of forcing and 
dissipation mechanisms [92) . This inertial range, with 



-^fein(^) ^ k~^/^, corresponds to merging and growth of 
vortices; the lEC and vortex growth can continue until 
the largest wavelength modes of the system contain sig- 
nificant energy in the form of a vortex (or vortex dipole) 
on the scale of the size of the system. The lEC is en- 
abled by the 2D nature of the system, where local vor- 
ticity must be normal to the 2D plane. This dimensional 
restriction gives rise to a conservation law for enstrophy 
(net squared vorticity) not present in 3D turbulence. As 
conserved energy and enstrophy flux cannot simultane- 
ously occur in two dimensions, 2D turbulence can simul- 
taneously display a flux of energy towards length scales 
larger than that of forcing (the lEC), and flux of en- 
strophy towards smaller length scales (a direct enstrophy 
cascade). The enstrophy cascade exhibits El-^{k) ex k~^ 
scaling over wavenumbers larger than that of the forcing 
scale. 

Whereas the classical lEC of a 2D fluid is associated 
with the growth of patches of vorticity due to vortex 
merging and pushing energy towards larger length scales, 
the enstrophy cascade is associated with the stretching of 
patches of vorticity. Experiments on forced and decaying 
2D turbulence have shown that an lEC that corresponds 
to vortex merging can appear with an enstrophy cascade 
[93] ■ or may be observed without a simultaneous enstro- 
phy cascade [Ml [SS] . Similarly, evidence for enstrophy 
cascades with approximate k~^ dependence in the en- 
ergy spectrum, and corresponding observations of vortex 
stretching, may appear without an lEC particularly in 
decaying 2D turbulence [Ml [HZ]- The type of forcing used 
in experiments may thus have a significant effect on the 
characteristics of observed spectra and flow dynamics. 

These cascades are only beginning to be explored in 2D 
QT. Among the most significant research topics within 
this field are the determination of (i) conditions under 
which an lEC and an enstrophy cascade can be found 
in forced or decaying 2D QT, either together or indi- 
vidually; (ii) superfluid dynamics that accompany either 
cascade process, in particular the clustering of quantized 
vortices accompanying an lEC, which is the process most 
readily envisioned to correspond to vortex growth in 2D 
QT, and the appearance of large-scale flow that may be 
produced by energy flux into large scales; and (iii) vortex 
dynamics that may provide insight into the mechanisms 
underlying turbulent dynamics, such as vortcx-antivortex 
annihilation and few-vortex dynamics. The remainder of 
this section is devoted to addressing aspects of these open 
problems. 

In 2D QT enstrophy is proportional to the number 
of vortices [551 ISE]- In decaying 2D QT, enstrophy is 
rarely conserved throughout the entire system, as vortex- 
antivortex annihilation events occur [981 I99j and serve 
as a dissipation mechanism. However, annihilation does 
not automatically preclude the appearance of an lEC in 
a forced system, and evidence for lECs have been ob- 
served in numerical simulations. These will be discussed 



below. For decaying 2D QT, existence of lECs remains 
an open question. To complicate matters further, in 
smaller systems, vortex dynamics are strongly influenced 
by the condensate confining geometry. For harmonically 
trapped BECs, the energy of a single vortex can be lost 
to the thermal cloud as the vortex precesses on a radial 
trajectory from the condensate centre towards the BEC 
boundary [45lfT00]. 

Another effect potentially prohibitive to establishing 
lECs in BECs is the generation of large amounts of 
phononic excitations. Sound can be generated from 
vortex-dipole recombination events, the movement of 
vortices, and in particular by the method applied to 
induce vortices into the condensates. Sound may also 
increase the frequency of vortex-antivortex annihilation 
events. Furthermore, in consideration of Onsager's argu- 
ments [1] , there may be a minimum local number density 
of vortices or a forcing rate necessary to generate an lEC 
process. These issues imply the presence of lEC processes 
in atomic BECs is no more clear cut than it is in classical 
turbulence, and considerable theoretical and experimen- 
tal studies have been concentrated on determining con- 
ditions for its existence. It is furthermore necessary to 
establish the dominant parameters that dictate different 
regimes of vortex dynamics for continuously forced and 
decaying 2D QT, so that these various regimes can be 
associated with features of energy spectra, forcing, and 
dissipation. 

To date, most of the work on compressible 2D QT has 
been numerical in nature. We turn our attention to sum- 
marizing some of the findings that are among the most 
relevant for understanding vortex dynamics and energy 
cascades. Decaying homogeneous 2D QT generated by 
random phase initial conditions was observed to exhibit 
a direct energy cascade with fc~^/^ scaling [55]. In a 
study of vortex clustering, statistical measures of clus- 
tering were applied to quantify the distribution of vor- 
tices nucleated from a moving laser stirrer in a harmon- 
ically trapped condensate. This study found no increase 
in clustering of vortices and only a significant degree of 
clustering was observed when the clustering was forced 
by the path of the laser stirrer 42 . In a different study 
using an annular trap, signs of weak clustering were ob- 
served, along with energy spectra consistent with (but 
not directly confirming) an lEC; this study was made in 
conjunction with experimental observations, and will be 
discussed at the end of this section [22j . Reeves et al [13] 
looked at a laser stirrer nucleating turbulent vortices and 
exciting acoustic energy into the BEC in order to distin- 
guish regimes of vortex and acoustic turbulence. It was 
found that in a particular regime of stirring where clus- 
tering of vortices of the same sign occurred, Kolmogorov 
k~^'^ scaling was observed, but became more transient 
as the fraction of vortices that were clustered decreased. 

Another recent 2D QT study of a homogeneous BEC 
simulated superfluid flow past past four static stirrers. 



and varied the degree of damping in the system. Regimes 
of vortex clustering and energy spectra were then deter- 
mined. For one range of damping, incompressible kinetic 
energy spectra approximately proportional to k~^'^ for 
length scales larger than that of the forcing scale were 
found |101| . Along with this spectral power law, a nega- 
tive flux of incompressible kinetic energy into the longest 
wavelength modes of the system was measured. The 
growth of vortex clusters was simultaneously observed: 
localized regions of vortices of the same sign were ob- 
served to grow in circulation quanta, and in size by a 
factor of approximately 5, both in time and distance away 
from the forcing region. Together, these three observa- 
tions show that lECs can be supported in 2D QT and 
that there is a correspondence between an lEC and the 
growth of patches of vorticity coinciding with vortex clus- 
ters, as in 2D classical turbulence. Nevertheless, there 
still remain numerous open issues regarding vortex clus- 
tering and lECs that must be resolved. These include 
a broad discovery of forcing and dissipation parameters 
that lead to the development of an TEC and vortex clus- 
tering, the roles and effects of boundary conditions and 
system geometry, the possibility of an enstrophy cascade 
developing, and the conditions under which a large region 
of vorticity on the scale of the system size might develop. 

There are indeed promising results towards observing 
a large vorticity end state in decaying 2D QT in recent 
numerical simulations of Simula et al. |102| . In this pre- 
liminary work, vortices were initially induced via phase 
imprinting, and the loss of enstrophy at system bound- 
aries was limited by a weak repulsive gaussian potential 
within the trap in addition to hard-wall boundary con- 
ditions for the trap. Clear separation and aggregation 
of vortices of like sign was observed on the scale of the 
condensate size. For the form of the trapping potential 
applied in this study, the minimum energy for a vortex 
within the condensate is at the condensate centre, in con- 
trast to harmonically trapped condensates, where it is 
a maximum. As in the previously discussed work, this 
demonstrates the significance of system parameters on 
observed phenomena. 

Finally, an analytical approach to understanding the 
spectra of 2D QT for a compressible superfluid was pur- 
sued in |28| . In this work, the shape of vortex cores in 
a BEC was used to analytically derive an expression for 
incompressible kinetic energy spectra of various vortex 
conflgurations in a homogeneous BEC. Intriguingly, for 
small enough length scales, the spectrum is proportional 
to fc^"^, although does this not correspond to an enstro- 
phy cascade. A particular type of vortex clustering was 
shown to lead to a k^^^^ distribution in the spectrum 
for large enough length scales. This paper proposed a 
scenario of forced 2D QT in which an lEC might appear 
over length scales larger than an assumed forcing scale 
at the lower end of the k~^ range, although there is no 
enstrophy cascade in the proposed picture. By consider- 



ing the fraction of vortices involved in clusters, numerical 
results of forced turbulence correspond surprisingly well 
to the analytical predictions |101j . An analytical deriva- 
tion of the Kolmogorov constant for forced 2D QT in a 
BEC was also proposed, suggesting that new approaches 
to 2D QT could eventually lead to new insights on 2D 
turbulence. 

From the body of recent theoretical work performed 
to date, we can conclude that the lEC, vortex clustering, 
and system-scale vortex growth may appear in 2D QT, 
but understanding the full range of conditions for the 
appearance of these features will require much further 
exploration. Simulations have shown that the manner in 
which vortices, as well as sound, are forced into the BEC 
can play a large role in the spectra observed, and this 
is one of the more complex aspects of the 2D QT field. 
In the above summaries, we have focused on investiga- 
tions that center on lECs and inertial regimes in energy 
spectra, but there are indeed other important aspects 
of 2D QT in BECs. In particular, discussions of non- 
thermal fixed points lead to new means of understanding 
and discussing quantum turbulence |25H27| . and studies 
of vortex dynamics and correlations may provide further 
means of analyzing regimes of QT ^\ . 



Experiments on two-dimensional QT 

As with 3D QT in BECs, there are still few exper- 
iments on 2D QT in trapped atomic condensates, and 
there is much to learn. The first 2D QT experiment 
by Neely et al.[22| observed the formation and evolu- 
tion of disordered distributions of vortices forced into a 
highly oblate annular BEC by stirring with a laser beam, 
as shown in Fig. [2jB). The technique excited very little 
acoustic energy, and injected on the order of 20 vortices, 
which may well be enough to display the physics of tur- 
bulence in such a small system. Subsequent to stirring, 
the system evolved into a persistent current as vortex 
number decayed; while this observation corresponds to 
the growth of energy in large-scale flow given small-scale 
forcing, the extent to which this observation directly re- 
lates to the vortex and spectral dynamics described above 
is currently unclear. Numerical simulations of the pro- 
cess showed the presence of pairs of like-signed vortices, 
an incompressible energy spectrum proportional to k~^'^ 
over a range of length scales larger than the forcing scale, 
and conservation of enstrophy over the time period as- 
sociated with the growth of the k~^/^ spectrum. While 
these observations are consistent with an lEC, the energy 
flux has not been determined, and the definite existence 
of an lEC in these simulations remains an open issue. 

Other methods have also been experimentally explored 
for producing 2D QT in highly oblate BECs [^, as 
shown in Fig. [2lC)-(F). These include modulating the 
strength of the trapping potential, suddenly applying 
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FIG. 2: (A). 2D QT can be studied in a highly oblate BEC, as 
shown in an in-situ image of the BEC taken along the tightly 
confining axis and a corresponding in-situ image taken along 
the radial plane (inset). An annular BEC (not shown) can be 
created by directing a repulsive focused laser beam through 
the BEC along the tight trapping direction. (B) In an an- 
nular BEC excited by small-scale stirring using a repulsive 
laser potential, vortices were observed after the repulsive cen- 
tral laser barrier was removed prior to BEC expansion. Note 
that the circulation of vortices can not be determined in this 
expansion method. Taken from the experiment of [21]. The 
remaining images of expanded BECs show 2D QT generation 
due to (C) modulation of the harmonic potential in a trap 
with a central focused repulsive laser potential, followed by 
50 ms of hold time, (D) laser beam modulation (image taken 
100 ms after removing beam), (E) sudden localized repulsive 
laser beam application and removal followed by 125 ms of hold 
time, and (F) trap spinning followed by 2.5 s of hold time. In 
all cases, the BEC is expanded prior to imaging. See Ref. [23| 
for details on these images and experimental methods. 



then removing a repulsive laser potential to a localized 
region in the BEC, modulating the intensity of a local- 
ized repulsive laser potential, and spinning a slightly el- 
liptical highly oblate trapping potential within the radial 
plane. Within the parameter ranges explored, all of these 
methods have been experimentally found to induce large 
numbers of vortices in harmonic and annular traps, and 
are candidates for further 2D QT studies. Unfortunately 
for the study of vortex turbulence, these methods also 
generally excite acoustic or breathing modes of the BEC, 
and may render 2D QT studies difficult without further 
modification. A deeper understanding the physical ori- 
gins of vortex generation in each case is also needed. 
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QUANTUM TURBULENCE IN BECS: FUTURE 
PROSPECTS 



Central experimental challenges 

Here we outline some of the main experimental chal- 
lenges that must be overcome to explore aspects of QT 
in BECs. Ideally, an experiment would be able to watch 
all vortices in real time, examine their dynamics and how 
they interact with the other vortices of a BEC, determine 
the circulations of all vortices, and measure correspond- 
ing energy spectra. This is an exceptionally challenging 
task! But there are likely to be realistic methods for 
approaching at least some of these goals. In 2D QT, 
methods for in situ imaging of vortices in the plane of 
a highly oblate BEC are currently being explored, and 
these show promise for use in multiple imaging methods 
that will allow for real-time probes of 2D QT vortex dy- 
namics. The 3D case is much more demanding, and will 
be briefly discussed in the next section. 

Atom interferometry presents an available method for 
circulation measurement [53J 11031 1104] . and in princi- 
ple may be used even in three dimensions. In practice, 
whether in two or three dimensions, the basic application 
of atom interferometry in these cases would involve the 
interference of a turbulent BEC with a reference BEC. 
Trapping two BECs simultaneously is challenging; driv- 
ing one and not the other into a turbulent state adds 
another layer of challenge that may be achievable with 
laser stirring beams, for example. Matter wave interfer- 
ence that is suitable for resolving circulations of a high- 
vorticity BEC is another open experimental challenge. 
Add to this the vortex tangles of 3D turbulence, and the 
experimental challenges are formidable, but not out of 
the question. 

Another experimental goal involves the generation of 
QT with minimal acoustic excitation so that forced QT 
can be reliably studied. Forcing techniques are currently 
one of the main topics of numerical investigation. As 
mentioned, a major experimental advance will occur with 
techniques that permit measurement of the kinetic en- 
ergy spectrum of a BEC, although the correspondence of 
the total kinetic energy spectrum with the incompress- 
ible component is an open question. Finally, on-demand 
vortex generation and manipulation techniques will al- 
low for studies of vortex dynamics and interactions, as 
mentioned earlier. 



ments, where solid Hydrogen tracer particles were ap- 
plied to visualize vortices in superfiuid ''He [52l 11071 fl08j 
the Monash University group of K. Helmerson is setting 
up an experiment in a two-component Bose-Einstein con- 
densate [102] . One component will be applied as tracer 
particles that sit inside vortices in the second component. 
The particles of the tracer component are imaged to track 
the vortex core dynamics in the turbulent condensate. 
This experiment aims to achieve real-time in-situ imag- 
ing of vortex dynamics and accurate determination of 
vortex core positions. Using this information it will be 
possible to reconstruct various aspects of the superfiuid 
flow. 



A recent theoretical proposal has shown impurity 
atoms can be applied to out-perform state of the art ther- 
mometry of condensates, and accurately measure tem- 
peratures in the nano-Kelvin regime |109) . Applying sim- 
ilar ideas, we propose a few impurity atoms sitting in a 
condensate could be used to directly measure the veloc- 
ity correlations of a turbulent vortex tangle |110j . Re- 
carti et al |105| showed that a single impurity atom can 
be coupled to a superfiuid reservoir via a Raman transi- 
tion. In a simple set-up, we consider two or more impu- 
rity atoms that are well separated and therefore do not 
interact with each other. When the coupling with the 
condensate is switched on, these probe impurity atoms 
pick up the classical correlations in the velocity field of 
the quantum vortex tangle. After turning off the cou- 
pling to the condensate, full quantum state tomography 
of the probe impurity atoms (as has been demonstrated 
for superconducting qubits |111U112| ) enables one to fully 
reconstruct the density matrix and therefore also directly 
reconstruct the velocity correlations. Such a setup allows 
in-situ, non-destructive measurements of properties of a 
turbulent vortex tangle. As these measurement could be 
potentially performed with a few well-separated probe 
atoms and the interaction with the condensate can be 
finely controlled, the dynamics of the vortex tangle will 
not be infiuenced by the measurement of the velocity cor- 
relations. Such direct measurements of turbulent velocity 
correlations should be distinguished from a reconstruc- 
tion of the velocity field, a posteriori, from knowledge of 
the vortex positions and will also detect any phononic 
correlations. 



Potential experimental probes of turbulence 



Impurity atoms are a potentially useful probe of Bose- 



Einstein condensates |105| I106J and a promising tool for 
visualizing and determining properties of quantum vor- 
tex tangles. Motivated by the approach used in the 
University of Maryland counterfiow turbulence experi- 



A second experimental approach to measure velocity 
correlations has been demonstrated in a different context, 
toward observation of an analog of the dynamical casimir 
effect for phononic excitations in a BEC of metastable he- 
lium atoms |113) . After trap modulation, the condensate 
is released and atom velocities are measured by a single 
atom detector that is also sensitive to the arrival position. 
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CONCLUSIONS 

To summarise, atomic BECs are a highly tuneable sys- 
tem that hold much promise for the development of the- 
oretical and experimental insights into some of the unan- 
swered questions surrounding the theory of quantum tur- 
bulence. This holds particularly true regarding aspects 
of compressibility and QT, and the nature of 2D QT, 
both of which have not been explored prior to recent 
work with BECs. With BECs, emerging methods for 
controlled vortex creation may soon allow deterministic 
preparation of initial states necessary for investigations 
into the chaotic dynamics of few-vortex systems. Rou- 
tine creation of highly oblate condensates provides the 
first system in which 2D quantum turbulence can be ex- 
perimentally explored. Atomic gas superfluids also pro- 
vide numerous other opportunities not mentioned in this 
paper, such as possibilities to investigate QT in spinor 
systems, or in degenerate Fermi gases. 

Atomic BECs are also currently the most accessible 
system in which to study the small scale properties of 
turbulent vortex flow in three dimensions. Theoretical 
investigations have begun to build up a picture of vortex 
dynamics and the processes contributing to the forcing 
and decay of turbulence at small scales and there is great 
scope for experimental verification in atomic BECs. At 
large scales, features of classical turbulence is an emer- 
gent feature of quantum turbulence, suggesting that re- 
search with BEC systems may provide insight into some 
of the outstanding questions of turbulence. 

BPA acknowledges funding from the US National Sci- 
ence Foundation, grant PHY-1205713. ACW acknowl- 
edges funding from EPSRC grant No. EP/H027777/1. 
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